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Abstract. We investigated the classical Shastry-Sutherland lattice under an external magnetic 
field in order to understand the recently discovered magnetization plateaux in the rare-earth 
tetraborides compounds RB4. A detailed study of the role of thermal fluctuations was carried 
out by mean of classical spin waves theory and Monte-Carlo simulations. Magnetization quasi- 
plateaux were observed at 1/3 of the saturation magnetization at non zero temperature. We 
showed that the existence of these quasi-plateaux is due to an entropic selection of a particular 
coUinear state. We also obtained a phase diagram that shows the domains of existence of 
different spin configurations in the magnetic field versus temperature plane. 



The Shastry-Sutherland lattice (SSL) was first considered in the early eighties by Shastry and 
Sutherland [Ij. This frustrated lattice can be described as a square lattice with antiferromagnetic 
couplings J' and additional diagonal antiferromagnetic couplings J in one square over two (figure 
d]). The quantum version of the SSL was widely studied both theoretically and experimentally 
after its experimental realization in the SrCu2(B03)2 compound. The quantum SSL was found 
to present a rich phase diagram as the ration J' /J varies. The existence of magnetization 
plateaux was theoretically predicted and experimentally verified for various rational values of 
the saturated magnetization M/Msat = 1/8, 1/4, 1/3 and 1/2 0. 

Recently magnetization plateaux were discovered in rare-earth tetraborides RB4 in which 
the rare-earth ions are placed in the (001) plane according to a lattice that is topologically 
equivalent to the SSL (figure [2]). RB4 compounds present large total angular momenta that 
justify a classical description of the SSL. High magnetic-fields experiments revealed plateaux in 
the magnetization curves of TbB4 at M/Msat = 1/4, 1/3 and 1/2 [3], in ErB4 at M/Msat = 1/2 
in and in TmB4 M/Msat = 1/2 and M/M^at = 1/8 [5]. 

We study Heisenberg spins on the SSL in the classical limit. The spins are represented by vectors 
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Figure 1. The SSL with 
magnetic couplings J' bonds 
along the edges of the squares 
and J along the diagonals. 



Figure 2. The topolog- 
ically identical structure 
realized in SrCu2 (603)2 
and in RB4. 



Figure 3. The VVD 

state at M/Msat = 1/3 
and the 12 spins unit cell 
(bold lines). Each triangle 
contains two spins up (red) 
and one spin down (blue). 



Si with ||Si|| = 1. The magnetic field is applied to the z— axis. The Hamiltonian is given by: 
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The sums are taken over the edges of the squares with antiferromagnetic couplings J' and over 
the diagonals with J couplings. In the absence of applied magnetic field the lowest energy 
configuration (LEG) is coplanar. It is Neel ordered for J'/J > 1 else it is a spiral state with 
an angle (/? = vr it arccos (^) between nearest-neighbour spins. Those two possible values for 
if create a supplementary discrete degeneracy above the continuous one P]. For the particular 
ratio of magnetic couplings J'/J = 1/2, the Hamiltonian Eq. ([T]) can be written as a sum on 
triangles up to a constant term: 
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Where Sa = J2ieA total magnetization of one triangle and Na is the number of 

triangles {Na = N). Minimizing the energy on a single triangle, one obtains the constraint: 



Sa = i^ with J'/J = 1/2 



(3) 



The LEG is obtained when this constraint is satisfied in every triangle. In the absence of 
magnetic field it is a coplanar 120° structure with a 12 spins unit cell. Saturation is reached 
at hsat = 9 J' while the magnetic field at M/Msat = 1/3 is /11/3 = 3 J' (with J'/J = 1/2). At 
/i = different spin configurations have the same classical energy: for example, the "umbrella 
configuration" in which the three sublattices always have the same projection on the z— axis for 
all magnetic field applied and the "Up-Up-Down" (UUD) state. In this latter configuration each 
triangle contains two spins up and one spin down (figure^]). The existence of a classical plateaux 
requires the LEG to be UUD [6J. Since energetic considerations do not favour the UUD state 



over the other configurations no magnetization plateaux are expected at zero temperature. 
At finite temperature the most favourable configuration is also determined by the density of 
states above the LEG. Under thermal fluctuations the classical spin wave spectrum of a particular 
state (or subset of states) can present soft modes whose energy in spin deviations is quartic in 
state of quadratic. Those modes contribute less to the free energy and the state with more soft 
modes is favoured. This phenomenon is known as Order by Disorder [7j. 

We computed the spin-wave spectrum of the UUD state. The coordinates of the 12 spins of the 
unit cell are rewritten 

Si(rO = (6f(r,),6r(r,),l-l/2((6f)2 + (6r)2)) (4) 

Up to second order in spin deviations applied on the 12 sublattices the Hamiltonian reads: 

n = Eo+ ^ V*(-k)AlV.(k) + ^W„ (5) 

]<.,v=x,y n>2 

Where is the classical energy and Hn ~ ©(e") and V*(-k) = (ef(-k), . . . , e5'2(-k)) . 
Eigenvalues of the matrix M should be determined by solving an order 12 polynomial which 
could not be done analytically. Its determinant is given by: 

Det{M) oc ( — 2 + 3 cos kx — 3 cos 2k,j: + cos Sk^ + cos ky)'^ (6) 

The expression for the determinant shows that lines of soft modes can be found. However since 
the dimension of this submanifold is only 1, in the thermodynamic limit, those lines of soft modes 
no longer play a role and they will not create a fall in the specific heat |8]. In a similar way 
as in the classical triangular lattice, in the SSL, the UUD state is selected by entropy selection 
without any full branch of soft modes, in contrast to the kagome case [9]. 

We performed MC simulations using the standard Metropolis algorithm on square samples of 
size L = ^/N = 6, 12, 18, 24 and 30 with periodic boundary conditions. System sizes are chosen 
so as to be commensurable with the 12 spins unit cell represented in figure O Magnetization, 
susceptibility x = dM/dh and specific heat were computed. The thermal equilibrium was 
reached by collecting data with up to ~ 10^ MC steps on an average number of five simulations 
per point. 

Quasi-plateaux are observed in the magnetization curves at M/Mgat = 1/3. Their width 
increases with the system size which confirms the fact that the its stabilization is not due 
only to the lines of soft modes. 

Figure H] shows the effect of the temperature on the quasi-plateaux. As the temperature 
decreases the quasi-plateaux tend to disappear which is in agreement with the analytical 
prediction. Their width increases with temperature until they are destroyed by temperatures 
higher than T ~ 0.1. The susceptibility presents two very sharp peaks around hij^ that become 
rounded and disappear as the temperature increases. The positions of those peaks were studied 
in order to determine qualitatively the domain on existence of the collinear phase in the (/i, T) 
plane (see below). We performed MC simulations on L = 12 systems in order to obtain the 
schematic phase diagram presented in figure [5l Numbers are attributed to the susceptibility 
peaks from low to high temperature in magnetic field scan (curves with circles) and from low to 
high magnetic field in temperature scan (curves with triangles). 

This phase diagram presents similarities with the ones of the classical the Kagome [9j and 
triangular lattices JOJ. In the high temperature region T > 0.1) the system enters the 
disordered paramagnetic phase. The low temperature region contains three phases. The nature 
of the phase transitions is still to be determined [llj. The collinear UUD phase is the LEC for 
values of the magnetic field below and until /11/3. As the temperature tends to zero this phase 




Figure 4. Susceptibility and magnetization Figure 5. Phase diagram of the classical SSL 
(inset) for a L = 18 system at J' / J = 1/2. at J' / J = 1/2. 

width converges to a single point at exactly hij^ = 3/2. The phases below and above UUD 
contain coplanar states for which quasi-long range order is expected. In the low field region each 
triangles contains one spins below the (x, y) plane while the remaining two spins raise until the 
UUD state. In high field region the down spin raises until saturation iUUU state). 
We also studied the effect of a small modification e of the magnetic coupling ratio J' / J = 1/2. 
Magnetization plateaux are also observed for J' / J = 1/2 it e. UUD is no longer the LEG at 
"very low temperature" (i.e. accessible with lowest-order thermal fluctuations). In state the 
LEG is a non-collinear umbrella. However at higher temperature UUD becomes the LEG and 
there still exists a quasi-plateau phase in the (/i, T) plane surrounded by the same not collinear 
phases as in the case J' / J = 1/2. 

In conclusion, we found magnetization quasi-plateaux at M/Msat = 1/3 in the classical SSL and 
showed that they are due to the entropic selection of the collinear UUD state. A phase diagram, 
that reminds the one of the classical triangular lattice, was established in the (/i, T) plane for the 
particular coupling ratio J' / J = 1/2. We also showed that, because of an entropic selection of 
the collinear "UUD" state, those quasi-plateaux survive under small variations e of the magnetic 
coupling around J' / J = 1/2 which disfavor energetically this collinear configuration. 
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